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Abstract 

We review recent theoretical advances in cold atom physics concentrating on strongly 
correlated cold atoms in optical lattices. We discuss recently developed quantum op- 
tical tools for manipulating atoms and show how they can be used to realize a wide 
range of many body Hamiltonians. Then we describe connections and differences to 
condensed matter physics and present applications in the fields of quantum comput- 
ing and quantum simulations. Finally we explain how defects and atomic quantum 
dots can be introduced in a controlled way in optical lattice systems. 
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1 Introduction 



Atomic physics experiments with quantum degenerate Bose and Fermi gases 
are characterized by the distinguishing features that we have (i) a detailed 
microscopic understanding of the Hamiltonian of the systems realized in the 
laboratory, and (ii) complete control of the system parameters via external 
fields. In particular, atoms can be trapped and their motion controlled in 
magnetic and optical traps, allowing, for example, the realization of quantum 
gases with different dimensionality at effectively zero temperature. In addition, 
atoms have many internal states which can be manipulated using laser light 
and can be employed as a probe of the gas properties, and their collisional 
properties can be tuned with magnetic and optical Feshbach resonances. 

In the early days of atomic BEC experiments [1,2,3,4,5], the main focus was to 
investigate condensate properties of matter waves like coherence, as described 
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theoretically by the Hartree-Fock-Bogoliubov mean field theory for weakly 
interacting quantum gases. More recently, emphasis has shifted to strongly 
interacting systems, which are much more in line with present interests in 
theoretical condensed matter physics. In particular, as first pointed out in 
Ref. [6], strongly interacting systems can be realized with cold atomic gases in 
optical lattices, i.e. periodic arrays of microtraps generated by standing wave 
laser light fields. This leads to Hubbard type lattice models, where atomic 
physics provides a whole toolbox to engineer various types of Hamiltonians for 
ID, 2D and 3D Bose and Fermi systems which can be controlled by varying 
external field parameters in a time dependent way. In addition, atomic physics 
provides systematic ways of loading these lattices with atoms. A prominent 
example is the Mott insulator-superfluid quantum phase transition with cold 
bosonic atoms, as first observed in the seminal experiment by I. Bloch and 
collaborators [7]. More generally, we expect that cold atoms in optical lattices 
will be developed in the coming years as a general quantum simulator of lattice 
models, allowing experimental insight into phase diagrams for certain classes 
of (toy) models (such as high-T c superconductivity) and for parameter regimes, 
where no rigorous theoretical approaches exist. In addition, new theoretical 
challenges appear in this context, for example, the study of time-dependent 
phenomena. Besides the condensed matter aspects, the engineered Hubbard 
models have direct application in quantum computing, where the controlled 
interactions can be used to create entanglement with high fidelity. 

In the present article we discuss the atomic physics point of developing a 
Hubbard toolbox by controlling laser and collisional interactions [8] . The ba- 
sic physical mechanisms for creating periodic optical trapping potentials are 
introduced in Sec. 2. In Sec. 3 we discuss some of the various Hubbard models 
that can be engineered in optical lattices and also give a microscopic explana- 
tion of the terms appearing in the Hamiltonians. Then we point out the simi- 
larities and in particular the novel aspects in comparison to condensed matter 
systems in Sec. 4, and proceed by describing possible applications in quantum 
computing and quantum simulations in Sec. 5. Finally, some prospects for in- 
troducing imperfections and quantum dots into optical lattices are discussed 
in Sec. 6 and we conclude in Sec. 7. 



2 Optical lattices 

Lasers are a very versatile tool for manipulating atoms. Atoms may be cooled 
or can even be trapped by laser light [9] . For trapping the purely conservative 
dipole force exerted by a laser with an inhomogeneous intensity profile is 
used [10]. In this section we present several examples of different laser-atom 
configurations which allow the realization of a variety of trapping potential, in 
particular periodic optical lattices with different geometries, and even trapping 
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potentials whose shape depends on the internal (hyperfine) state of the atom. 
In our derivations we will neglect spontaneous emission and later establish the 
consistency of this approximation by giving an estimate for the rate at which 
photons are spontaneously emitted in a typical optical lattice setup. 



2. 1 Optical potentials 



The Hamiltonian of an atom of mass m is given by Ha = p 2 /2m+J2j \ e j) ( e j I 
Here p is the center of mass momentum operator and \ej) denote the inter- 
nal atomic states with energies Uj (setting h = 1). We assume the atom to 
initially occupy a metastable internal state |e ) = \a) which defines the point 
of zero energy. The atom is subject to a classical laser field with electric field 
E(x, t) = E(xl, t)eexp(— iuot) where u is the frequency and e the polarization 
vector of the laser. The amplitude of the electric field -E(x, t) is varying slowly 
in time t compared to 1/w and slowly in space x compared to the size of the 
atom. In this situation the interaction between atom and laser is adequately 
described in dipole approximation by the Hamiltonian H^ p = — /iE(x, t)+h.c, 
where fi is the dipole operator of the atom. We assume the laser to be far 
detuned from any optical transition so that no significant population is trans- 
ferred from | a) to any of the other internal atomic states via i^dip- We can 
thus treat the additional atomic levels in perturbation theory and eliminate 
them from the dynamics. In doing so we find the AC-Stark shift of the internal 
state \a) in the form of a conservative potential V(x.) whose strength is deter- 
mined by the atomic dipole operator and the properties of the laser light at 
the center of mass position x of the atom. In particular V(x) is proportional 
to the laser intensity \E(k, t)\ 2 . Under these conditions the motion of the atom 
is governed by the Hamiltonian H = p 2 /2m + V(x). 

Let us now specialize the situation to the case where the dominant contribu- 
tion to the optical potential arises from one excited atomic level |e) only. In 
a frame rotating with the laser frequency the Hamiltonian of the atom is ap- 
proximately given by Ha = p 2 /2m + 8 |e) (e| where 5 = u e — u is the detuning 
of the laser from the atomic transition |e) <-> \a) . The dominant contribution 
to the atom-laser interaction neglecting all quickly oscillating terms (i.e. in the 
rotating wave approximation) is given by ifdip = ^(x) |e) (a\ /2 + h.c. Here 
Q = — 2£'(x, t) (e| fie \a) is the so called Rabi frequency driving the transitions 
between the two atomic levels. For large detuning 5 ^> Q adiabatically elim- 
inating the level |e) yields the explicit expression V(x) = |fi(x)| 2 /45 for the 
optical potential. The population transferred to the excited level |e) by the 
laser is given by |f2(x)| 2 /4<5 2 and this is the reason why we require f2(x) <C 5 
for our adiabatic elimination to be valid. 
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Fig. 1. Band structure of an optical lattice of the form Vq(x) = Vbcos 2 (/cx) for 
different depths of the potential, a) Vq = 5Er, b) Vq = IOEr, and c) Vb = 25Er. 

£.1? Periodic lattices 

For creating an optical lattice potential we start by superimposing two counter 
propagating running wave laser beams with £±(x, t) = exp(±zfoc) prop- 
agating in x-direction with amplitude E , wave number k and wave length 
A = 2ir/k. They create an optical potential V(~x) oc cos 2 (A;x) in one dimension 
with periodicity a = A/2. Using two further pairs of laser beams propagating 
in y and z direction respectively, a full three dimensional periodic trapping 
potential of the form 

y(x) = V 0x cos 2 (kx) + V 0y cos 2 (%) + V 0z cos 2 (kz) (1) 

is realized. The depth of this lattice in each direction is determined by the 
intensity of the corresponding pair of laser beams which is easily controlled in 
an experiment. 

Bloch bands and Wannier functions: For simplicity we only consider one spatial 
dimension in Eq. (1) and write down the Bloch functions (j)^(x) with q the 
quasi momentum and n the band index. The corresponding eigenenergies M n ' 
for different depths of the lattice V /E R in units of the recoil energy E R = 
are shown in Fig. 1. Already for a moderate lattice depth of a few 
recoil the separation between the lowest lying bands is much larger than their 
extend. In this case a good approximation for the gap between these bands is 
given by the oscillation frequency u>t of a particle trapped close to one of the 
minima Xj (= lattice site) of the optical potential. Approximating the lattice 
around a minimum by a harmonic oscillator we find lot = V^VoEr [6] . 

The dynamics of particles moving in the lowest lying well separated bands will 
be described using Wannier functions. These are complete sets of orthogonal 
normalized real mode functions for each band n. For properly chosen phases 
of the (j)^'{x) the Wannier functions optimally localized at lattice site 
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Fig. 2. a) Simple three dimensional cubic lattice, b) Sheets of a two dimensional 
triangular lattices, c) Set of one dimensional lattice tubes. In a)-c) tunnelling of 
atoms through the optical potential barriers is only possible between sites which 
are connected by lines. 

defined by [11] 

w n (x-x j ) = e-v 2 J2 e ~ iqXj <f> < ? ) (xl ( 2 ) 

where is a normalization constant. Note that for Vo — > oo and fixed k the 
Wannier function w n (x) tends towards the wave function of the n-th excited 
state of a harmonic oscillator with ground state size a = ^Jl/muT- We will use 
the Wannier functions to describe particles trapped in the lattice since they 
allow (i) to attribute a mean position Xj to the particles in a given mode and 
(ii) to easily account for local interactions between particles since the dominant 
contribution to the interaction energy arises from particles occupying the same 
lattice site Xj. 

Lattice geometry: The lattice site positions determine the lattice geometry. 
For instance the above arrangement of three pairs of orthogonal laser beams 
leads to a simple cubic lattice (shown in Fig. 2a). Since the laser setup is very 
versatile different lattice geometries can be achieved easily. As one example 
consider three laser beams propagating at angles 2tt/3 with respect to each 
other in the xy-plane and all of them being polarized in z direction. The 
resulting lattice potential is given by V(x) oc 3 + 4 cos(3fcx/2) cos( v / 3%/2) + 
2 cos(v / 3%) which is a triangular lattice in two dimensions. An additional pair 
of lasers in z direction can be used to create localized lattice sites (cf. Fig. 2b). 
Furthermore, as will be discussed later in Sec. 3, the motion of the atoms can 
be restricted to two or even one spatial dimension by large laser intensities. 
As will be shown it is thus possible to create truly one- and two-dimensional 
lattice models as indicated in Fig. 2b, c. 

Pseudo random site offset: This can be achieved by additional superlattice 
potentials. Two additional laser beams propagating at small angles a and at 
7r — a with respect to the x- axis of the lattice, respectively, superimpose a 
slowly varying potential with spatial period I = sin(a)/2A; where k is the wave 
number of the additional lasers. If I 3> a the main effect of the superlattice 
is an energy offset q oc cos(xi/l) for lattice site Xj. If the periodicity I is 
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Fig. 3. a) Atomic fine and hyperfine structure of the most commonly used alkali 
atoms 23 Na and 87 Rb. b) Schematic AC-stark shift of the atomic level S\/2 with 
m s = 1/2 (dashed curve) and with m s = — 1/2 (solid curve) due to the laser beam 
<7+ as a function of the laser frequency u>. The ac-stark shift of the level S\/2 with 
m s = —1/2 can be made by choosing the laser frequency u = ujl- 

incommensurate with a and several additional lasers with different periodicity 
l n are superimposed to the original lattice the site offsets q can be made quasi 
random over the size of the actual optical lattice. 

State dependent lattices: As already mentioned above the strength of the op- 
tical potential crucially depends on the atomic dipole moment between the 
internal states involved. Thus we can exploit selection rules for optical transi- 
tions to create differing traps for different internal states of the atom [12,13,14]. 
We will illustrate this by an example that is particularly relevant in what fol- 
lows. We consider an atom with the fine structure shown in Fig. 3a, like e.g. 
23 Na or 87 Rb, interacting with two circularly polarized laser beams. The right 
circularly polarized laser a + couples the level Si/ 2 with m s = — 1/2 to two ex- 
cited levels P1/2 and P 3 / 2 with m s = 1/2 and detunings of opposite sign. The 
respective optical potentials add up. The strength of the resulting AC-Stark 
shift is shown in Fig. 3b as a function of the laser frequency to. For u = ujl the 
two contribution cancel. The same can be achieved for the a~ laser acting on 
the 5*1/2 level with m s = 1/2. Therefore at lu = the AC-Stark shifts of the 
levels Si/2 with m s = ±1/2 are purely due to a± polarized light which we de- 
note by V±(x). The corresponding level shifts of the hyperfine states in the S1/2 
manifold (shown in Fig. 3a) are related to V±(x) by the Clebsch-Gordan co- 
efficients, e.g. V| F=2 , mF= 2) (x) = V+(x), V ]F=hmF=1) (x) = 3V + (ar)/4 + V_(x)/A, 
and V\ F =i, mF =-i) (x) = V+(x)/4 + 3V_(z)/4. 

State selectively moving the lattice: The two standing waves a± can be pro- 
duced out of two running counter-propagating waves with the same intensity 
as shown in Fig. 4. Moreover it is possible to move nodes of the resulting stand- 
ing waves by changing the angle of the polarization between the two running 
waves [12]. Let {ei,e 2 ,e 3 } be three unit vectors in space pointing along the 
{x, y, z] direction, respectively. The position dependent part of the electric 
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Fig. 4. Laser configuration for a state selective optical potential. Two standing circu- 
lar polarized standing waves are produced out of two counter-propagating running 
waves with an angle 2<p between their polarization axes. The lattice sites for different 
internal states (indicated by closed and open circles) are shifted by Ax = 2(p/k. 

field of the two running waves Ei j2 is given by Ei oc e tkx (cos(^)e3 + sin(<£>)e 2 ), 
E 2 oc e~ tkx (cos(<^)e3 — sin(<£>)e 2 ). The sum of the two electric fields is thus 
E x + E 2 oc cos(fcr — (p)cr_ — cos(kx + <p)cr + , where <j± = e 2 ± ie 3 and the 
resulting optical potentials are given by 



By changing the angle (p it is therefore possible to move the nodes of the two 
standing waves in opposite directions. Since these two standing waves act as 
internal state dependent potentials for the hyperfine states the optical lattice 
can be moved in opposite directions for different internal hyperfine states. 

Population transfer between hyperfine states: Two atomic hyperfine levels can 
be coupled coherently via magnetic dipole moment matrix elements by an 
additional oscillating microwave field. Alternatively transitions can also be 
driven by a Raman laser setup which consists of two laser beams and couples 
two hyperfine levels via a far detuned excited state. In both cases the dynamics 
is described by a Hamiltonian of the form H R = (Q R \a) (b\ +h.c.) /2 + 5 R \b) (b\ 
with fl R the Rabi frequency, 5 R the detuning and \a) , \b) the two coupled 
atomic states. 



2.3 Validity 

In all of the above calculations we have only considered the coherent interac- 
tions of an atom with laser light. Any incoherent scattering processes which 
lead to spontaneous emission were neglected. We will now establish the va- 
lidity of this approximation by estimating the mean rate T e g of spontaneous 
photon emission from an atom trapped in the lowest vibrational state of the 
optical lattice. This rate of spontaneous emission is given by the product of 
the life time T of the excited state and the probability of the atom occupying 
this state. In the case of a blue detuned optical lattice (dark optical lattice) 
5 < the potential minima coincide with the points of no light intensity and 
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we find the effective spontaneous emission rate T e g m If the lattice 

is red detuned (bright optical lattice) 5 > the potential minima match the 
points of maximum light intensity and we find r cff ps TVq/5. Since Vo > 
the spontaneous emission in a red detuned optical lattice will always be more 
significant than in a blue lattice, however, as long as Vo <^ 5 spontaneous 
emission does not play a significant role. 

Typical numerical values: In a typical blue detuned optical lattice with A = 
514nm for 23 Na atoms (5*1/2 — -P3/2 transition at A 2 = 589nm and T = 2n x 
10MHz) the recoil energy is E R ps 2ir x 33kHz and the detuning from the 
atomic resonance is 5 ~ —2.3 x 10 9 -Er. A lattice with depth Vq = 25Er leads 
to a trapping frequency of ujt = 1-0 En yielding a spontaneous emission rate 
of r e fr ~ 10~ 2 /s while experiments are typically carried out in times shorter 
than one second. Therefore spontaneous emission does not play a significant 
role in such experiments. 



3 The (Bose) Hubbard model 



We consider a gas of interacting particles moving in an optical lattice. Starting 
from the full many body Hamiltonian including local two particle interactions 
we first give a naive derivation of the Bose-Hubbard model (BHM) [15,6] and 
present related models which can be realized in an optical lattice. Then we 
proceed by discussing adiabatic and irreverisble schemes for loading the lat- 
tice with ultracold atoms. Finally we examine the microscopic origin of the 
interaction terms appearing in the BHM. Throughout we will mostly concen- 
trate on bosonic atoms, however, similar derivations can also be performed for 
fermions [14,16] and we will also give one example of such a fermionic model. 



3. 1 Naive derivation of the BHM 



The Hamiltonian of a weakly interacting gas in an optical lattice is 

#fuii = / (Fx&(x) (j^ + V (x) + V T (^ *(x)+| J <fe#t(x)£t( x )£( x )£( x ) 

( 4 ) 

with \l> (x) the bosonic field operator for atoms in a given internal atomic 
state I b) and Vr(x) a (slowly varying compared to the optical lattice Vo(x)) 
external trapping potential, e.g. a magnetic trap or a superlattice potential. 
The parameter g is the interaction strength between two atomic particles. If 
the atoms interact via s-wave scattering only it is given by g — Ana s /m with a s 
the s-wave scattering length. We assume all particles to be in the lowest band 



8 



of the optical lattice and expand the field operator in terms of the Wannier 
functions *&(x) = J2i biW^(x — Xj), where b\ is the destruction operator for a 
particle in site X;. We find H fnti = - Ejj Jijb\bj + \ T,i,j,k,i ^ ijkhW-b k b u where 

Jij = ~ J dxw (x - Xj) + Vo(x) + Vr(x)^j u? (x - X,), 



and 



^ijfci = 5'/ dxw (x - Xj)w (x - x j )w (x - x fe )w (x - xj). 



The numerical values for the offsite interaction matrix elements Uijki involving 
Wannier functions centered at different lattice sites as well as tunnelling matrix 
elements J^ to sites other than nearest neighbors (note that diagonal tunneling 
is not allowed in a cubic lattice since the Wannier functions are orthogonal) 
are small compared to onsite interactions t/ ooo = U an d nearest neighbor 
tunneling J 01 = J for reasonably deep lattices V > 5E R . We can therefore 
neglect them and for an isotropic cubic optical lattice arrive at the standard 
Bose-Hubbard Hamiltonian 

#bh = - J E + \ E Whh + E <,'>?,- (5) 

Here (i, j) denotes the sum over nearest neighbors and the terms 6j = Vr(xj) 
arise from the additional trapping potential. The physics described by Hbh is 
schematically shown in Fig. 5a. Particles gain an energy of J by hopping from 
one site to the next while two particles occupying the same lattice site provide 
an interaction energy U. An increase in the lattice depth Vq leads to higher 
barriers between the lattice sites decreasing the hopping energy J as shown in 
Fig. 5b. At the same time two particles occupying the same lattice site become 
more compressed which increases their repulsive energy U (cf. Fig. 5b). 

Tunneling term J: In the case of an ideal gas where U = the eigenstates of 
i^BH are easily found for 6j = and periodic boundary conditions. From the 
eigenvalue equation = — 2J cos(qa) we find that 4 J is the height of the 
lowest Bloch band. Furthermore we see that the energy is minimized for q = 
and therefore particles in the ground state are delocalized over the whole lat- 
tice, i.e. the ground state of N particles in the lattice is \^sf) oc (J2i b\) N |vac) 
with | vac) the vacuum state. In this limit the system is superfluid (SF) and 
possesses first order long range off diagonal correlations [15,6]. 

Onsite interaction U : In the opposite limit where the interaction U dominates 
the hopping term J the situation changes completely. As discussed in detail 
in [15,6] a quantum phase transition takes place at about U ~ 5.8zJ where z 
is the number of nearest neighbors of each lattice site. The long range correla- 
tions cease to exist in the ground state and instead the system becomes Mott 
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Fig. 5. a) Interpretation of the BHM in an optical lattice as discussed in the text, 
b) Plot of scaled onsite interaction U/Er multiplied by a/a s (S> 1) (solid line with 
axis on left-hand side of graph) and J/Er (dashed line, with axis on right-hand side 
of graph) as a function of Vq /Er = V X} y jZ o / Er (for a cubic 3D lattice) . 

insulating (MI). For commensurate filling of one particle per lattice site this 
MI state can be written as \^mi) oc Jlj bj |vac) . 

3.2 Lattice geometry and spin models 

The lattice geometry and dimensionality are reflected in the sum over the 
nearest neighbors Exotic geometries like e.g. a Kagome lattice can be 

realized [17] and disordered systems may be studied by superimposing a pseu- 
dorandom potential [18]. If for the simple cubic lattice the laser intensities in z 
and y direction are made very large compared to the intensity in x tunnelling 
predominantly happens along the x direction (cf. Fig. 5) and the system is 
effectively reduced to a set of one dimensional lattices. In addition, when the 
laser intensities are turned on sufficiently slowly (adiabatically), the correla- 
tions between atoms in different optical lattice tubes cease via the SF to MI 
transition and the experimental setup becomes an ensemble of independent one 
dimensional systems. These can be used to study properties of strongly corre- 
lated gases like e.g. a To nks- Girardeau gas recently experimentally achieved 
in [19], or the excitation spectrum of one dimensional ultracold gases [20]. 

Spin models: If several species of atoms, e.g. differernt hyperfine states, are 
trapped in the lattice one often identifies them with spin degrees of freedom. 
The most common example is a spin 1/2 system created from two states 
||) = | a) and ||) = \b) . Such systems possess rich phase diagrams [21] and 
as shown in [22,23,24] one can engineer spin exchange interactions and design 
their properties like anisotropy and sign by proper choices of optical potentials. 
It is even possilbe to create Cooper pairs of bosonic atoms and realize the 
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regime of a Luttinger liquid [25]. 

For a lattice of fermionic atoms, using similar setups as discussed above and 
an additional Raman laser fl R , one can realize spin dependent Fermi Hubbard 
models [16,14] with e.g. a Hamiltonian 



where c a j (c a j) destroys (creates) a fermion in lattice site j and internal state 
a £ {t> !} an d the hopping matrix element J a can be made spin dependent 
by state selective optical potentials. Also loading of mixtures of bosonic and 
fermionic atoms into an optical lattice was proposed, the phase diagram was 
worked out [26] and effects like pairing of fermions with bosons [27] or the 
phase separation of the two species [28] were studied. 

The Hofstaedter butterfly: A slightly more involved laser setup even allows for 
the realization of effective magnetic fields in two dimensional optical lattices 
[29]. This requires time reversal symmetry to be broken which can e.g. be 
achieved by a combination of accelerating a state dependent lattice and using 
Raman lasers to induce the hopping between sites as schematically shown 
in Fig. 6. A particle moving around one plaquette of the lattice acquires a 
phase shift that is purely determined by the laser phases and thus can be 
chosen freely between 2na G [0, 2n] [29]. It corresponds to a magnetic flux of 
$ = 2ira j e (with e the electron charge) and the dynamics is therefore identical 
to an electron moving in a lattice subject to an external magnetic field. Any 
strength of the magnetic field - an in particular also very large field strengths 
- can be achieved in the optical lattice setup and thus it should be possible to 
obtain the fractal energy band structure (the Hofstaedter butterfly) predicted 
by Hofstaedter for lattice electrons moving in huge magnetic fields [30]. Such 
optical lattice setups also allow for investigations of fractional quantum Hall 
states [31]. 

In this article we exclusively discuss the case of tight optical confinement in all 
three dimensions with few particles per lattice site and the resulting Hubbard 
models. The situation where only one or two pairs of laser beams are used 
to create arrays of coupled BECs is instead similar to arrays of Josephson 
junctions. Such setups have been investigated intensively both theoretically 
[32,33,34,4] as well as experimentally [35,36]. Also early experiments which 
used purely optical cooling and loading methods to prepare atoms in the low- 
est Bloch band of the lattice and hence achieved small filling factors of approx- 
imately 10% are not considered here. For such fillings interaction effects are 
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Fig. 6. Optical lattice setup: Open (closed) circles denote atoms in state \a) (\b)). 
a) Hopping in the y-direction is due to kinetic energy and described by the hop- 
ping matrix element J v being the same for particles in states \a) and Along 
the x-direction hopping amplitudes are due to the additional lasers f^i and b) 
Trapping potential in x-direction. Adjacent sites are set off by an energy A because 
of the acceleration or a static inhomogeneous electric field. The laser fii is resonant 
for transitions | a) «-> \b) while O2 is resonant for transitions between \b) <-> \a) due to 
the offset of the lattice sites. Because of a spatial dependence of f2i,2 in y-direction 
atoms hopping around one plaquette get phase shifts of 2na. 

not important and one particle effects like Bloch oscillations, quantum chaos, 
and laser cooling (see e.g. [37,38]) were of main interest in these experiments. 



3.3 Loading schemes 



Only by using atomic Bose-Einstein condensates (BEC) has it become possi- 
ble to achieve large densities corresponding to a few particles per lattice site. 
A BEC can be loaded from a magnetic trap into a lattice by slowly turning on 
the lasers and superimposing the lattice potential over the trap. The system 
adiabatically undergoes the transition from the SF ground state of the BEC to 
the MI state for a deep optical lattice [6]. Note that the adiabaticity condition 
can easily be fulfilled in this scenario since the changing optical lattice po- 
tential does not induce significant quasi-momenta and therefore the relevant 
energy is the distance to the excited band. Experimentally this loading scheme 
and the SF to MI transition was realized in several experiments [7,20,39] and 
led to a MI with up to two particles per lattice site. The number of defects in 
the created 'optical crystals' were limited to approximately 10% of the sites. 

Defect suppressed optical lattices: One method to further decrease the number 
of defects in the lattice was recently described in [40]. The proposed setup 
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Fig. 7. a) Avoided crossings in the energy eigenvalues E for n = 1,2,3 (Note the 
variation in the vertical scale). For the chosen values of 5, and 5f (dotted vertical 
lines) only one avoided crossing is traversed transferring exactly one atom from \a) 
to | b) as schematically shown in b). c) Loading of an atom into the first Bloch band 
from a degenerate gas and spontaneous decay to the lowest band via emitting a 
phonon into the surrounding degenerate gas. 



utilizes two optical lattices for internal states \a) and \b) with substantially 
different interaction strengths Ub ^ U a and identical lattice site positions Xj. 
Initially atoms in \a) are adiabatically loaded into the lattice and brought 
into a MI state where the number of particles per site may vary between 
n = 1, ...n max because of defects. Then a Raman laser with a detuning varying 
slowly in time between 5i and 5/ is used to adiabatically transfer exactly 
one particle from \a) to \b) as shown in Fig. 7a,b. This is done by going 
through exactly one avoided crossing. During the whole of this process transfer 
of further atoms is blocked by interactions. This scheme allows a significant 
suppression of defects and - with additional site offsets q - can also be used 
for patterned loading [39] of the \b) lattice. 



Irreversible loading schemes: Further improvement could be achieved via ir- 
reversible loading schemes. An optical lattice is immersed in an ultracold de- 
generate gas from which atoms are transferred into the first Bloch band of 
the lattice. By spontaneous emission of a phonon [41] the atom then decays 
into the lowest Bloch band as schematically shown in Fig. 7. By atom-atom 
repulsion further atoms are blocked from being loaded into the lattice. This 
scheme can also be extended to cool atomic patterns in a lattice. In contrast 
to adiabatic loading schemes it has the advantage of being repeatable without 
removing atoms already stored in the lattice. 
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3.4 Microscopic picture 



To obtain a better understanding of the interaction term U we consider two 
atoms trapped in one lattice site neglecting the tunnelling term J. We ap- 
proximate the lattice potential by a harmonic oscillator which decouples the 
center of mass motion from the relative motion. While the center of mass co- 
ordinate moves in a harmonic potential of frequency uot the potential of the 
relative coordinate r is the sum of a harmonic trap with frequency ojt and 
the interatomic interaction potential V int (r) as schematically shown in Fig. 8. 
The potentials VJ nt (r) are well known and one can calculate the energy lev- 
els ei with good precision [42,43,44]. The interaction matrix element U can 
then be identified with the difference in energy U = E\ — 3out/2 of the level 
I which is closest to the ground state of the bare harmonic oscillator level. 
The states in the Bose-Hubbard picture are thus related to the microscopic 
picture as indicated in Fig. 8. Also inelastic processes leading to particle loss 
can be treated in detail and a microscopic understanding of the behavior of 
U close to magnetic and optical Feshbach resonances can also be gained from 
this picture [42,45,46]. Based on the detailed knowledge of the two particle 
levels quantum optical schemes for the photoassociation of two atoms in \a) 
to form a molecule \m) can be devised. For instance, as schematically shown 
in Fig. 8, using a set of Raman lasers for molecule creation from a two particle 
MI state was proposed in [47]. In comparison to customary photoassociation 
the motion of the atoms and the molecule are precisely controlled and the 
resulting dynamics is described by a Hamiltonian of the form 



with f2 m the effective Rabi frequency, 5 m the detuning and m {fh)) destroy- 
ing (creating) a molecule in state \m) . As above the mode \m) appearing 
in the many particle Hamiltonian if mo i has a well defined microscopic origin 
(cf. Fig. 8). In the case of three particles occupying one lattice site analytic 
understanding can be obtained from Efimov states which have been used to 
calculate three particle loss rates and interaction properties near Feshbach 
resonances [48,49,50]. 



4 Condensed matter aspects 

The setups discussed above realize - for the first time - strongly correlated 
systems of neutral atoms. They are described by Hamiltonians which have 
been subject to extensive investigation in condensed matter physics where 
numerous techniques were developed to find out about their properties. It 
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Fig. 8. Microscopic picture: The interatomic interaction plus optical potential V (not 
to scale) for the relative motion of two atoms. Levels above V = are identified 
with two atoms in the lattice and states below V = as molecular states since they 
differ significantly in their extend. The relation between levels in the Bose-Hubbard 
picture and in the microscopic picture are indicated by arrows. Note that we assume 
the optical lattice to be identical for atoms \a) and the molecular state \m) which 
is only true for highly excited molecular states. 

is not the purpose of this article to review these techniques but we rather 
point out the most important differences between customary condensed matter 
systems and atoms in optical lattices. 

Controllability: The theoretical models for atoms in lattices follow from a mi- 
croscopic picture, the origin of every term in the Hamiltonian is well justified 
and can be controlled via the experimental setup. For instance, the lattice ge- 
ometry and tunnelling barriers are controlled by lasers parameters and their 
arrangement. Additional flexibility comes from the ability to tune the interac- 
tion properties of the atoms with magnetic and optical Feshbach resonances 
[42,45,46]. An even larger class of condensed matter systems can be realized by 
exploiting the lattice as a universal quantum simulator as shown in Sec. 5.2. 

Inhomogeneity: The size of an optical lattice is relatively small in comparison 
to its periodicity a with typically a few hundred lattice sites in each dimension 
only. Also, the lattice is usually not enclosed in a box but it is put into a slowly 
varying harmonic potential which translates into a space dependent site offset 
(can equally well be viewed as a local chemical potential). This leads e.g. to 
soft boundaries of the system (cf. Fig. 9a), or the simultaneous coexistence of 
spatially separated phases as shown in Fig. 9b where alternating SF and MI 
phases are present. 

Unitary dynamics: One of the most striking differences is the excellent iso- 
lation of an optical lattice from its environment. As shown above for spon- 
taneous emission, one of the most dominant decoherence mechanisms, the 
typical coherence time is on the order of seconds. The characteristic ener- 
gies in the Hamiltonian lead to a unitary evolution on millisecond time scales 
and the external parameters like the laser intensity can be changed in even 
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shorter times. This allows the experimental study of the unitary dynamics of 
strongly correlated atoms from the adiabatic to the sudden regime [7,51], a 
topic usually not considered in traditional condensed matter physics. To ob- 
tain a better understanding of these dynamical processes new analytical and 
numerical methods have to be developed. Exact numerical calculations for 
small systems have been performed [47,52] and also trajectory based simula- 
tion methods were developed [53,54]. By combining techniques from quantum 
information theory with known properties of one dimensional systems [55] it 
was possible to develop a numerical method, the so called time evolving block 
decimation algortihm (TEBD) which is related to density matrix renormaliza- 
tion group (DMRG) techniques [56,57], for simulating the unitary dynamics of 
a large class of one dimensional lattice models. In Fig. 9c we show results from 
dynamical calculations on the SF to MI transition which indicate differences 
from static ground state calculation in the off-diagonal elements of the one 
particle density matrix. 

Quantum State engineering: Simple condensed matter Hamiltonians like the 
BHM are often considered as toy models for obtaining insight into more com- 
plicated real physical systems and are thus not expected to fully describe the 
physics and in particular the full many body wave function. Through their ex- 
act realization with neutral cold atoms and the excellent isolation of these sys- 
tems from their environment such toy models become interesting for quantum 
state engineering and quantum computing applications which require precise 
knowledge of the system wave function. Quantum optical methods provide 
the accurate control and manipulability and the strongly correlated nature 
of the system can e.g. be exploited to initialize a quantum register [6,12], 
achieve increased stability [58] or even to provide the basic resource for quan- 
tum computations [59]. Some of these quantum computing implementations 
will be presented in Sec. 5. 



5 Quantum computing implementations 

One of the big challenges for 21st century physics is the demonstration of a 
quantum computer which promises a significant speed up for certain kinds of 
algorithms, most prominently the factorization of large numbers and searching 
in unsorted databases [60]. While the realization of a universal quantum com- 
puter is certainly a long term goal, also a number of nontrivial applications 
with limited quantum computer resources exist for the nearer future in the 
areas of quantum communication, entanglement creation (see Sec. 5.1), special 
purpose quantum computing, or Feynman's universal quantum simulator (see 
Sec. 5.2). 

The basic requirements for a physical system to be capable of quantum com- 



16 




Fig. 9. Absolute values of the one particle density matrix p m n = (b m b n ) as a function 
of site indices m, n for a one dimensional lattice superimposed by an additional 
slowly varying harmonic oscillator potential. All results are obtained by the TEBD 
[51]. (a) SF ground state for U/2J = 2, (b) Intermediate ground state for U/2J = 6. 
This state clearly indicates the coexistence of MI and SF states in different regions 
of the system. The SF region in the center of the trap is surrounded by two MI 
regions with one particle per lattice site followed by a small superfluid at the edges 
of the trap, (c) Intermediate state created from the SF by dynamically ramping up 
the depth of the lattice until U/2J = 6 [51]. Note the correlations between the two 
SF regions which are absent in the corresponding ground state. 



puting have been laid out by DiVincenzo [61]. Many of these are fulfilled for 
atoms in optical lattices. An optical lattice provides an intrinsically scalable 
and well defined set of qubits (quantum register) with long decoherence times. 
These are usually constituted by two metastable ground states of the atoms \a) 
and | b) trapped in the lattice. The initialization of the quantum register can 
be achieved by the loading schemes discussed above where the MI state with 
one particle per lattice site corresponds to the register initialized with zeros. 
Any quantum computation can be carried out by applying single qubit gates 
via Raman lasers, and controlled interactions between the particles to per- 
form two qubit gates as described in Sec. 5.1 [12,13,62,63]. Also several other 
schemes for quantum computing in optical lattices have been proposed based 
on e.g. strong dipole-dipole interactions between Rydberg atoms [64], the mo- 
tional state of the atoms [65], atom tunnelling between neighboring sites [66], 
or measurements on highly entangled states [59]. Currently the main (techni- 
cal) drawback in many of these optical lattice quantum computing schemes is 
the difficulty of addressing individual atoms. While a large number of quan- 
tum gate operations can already be carried out in parallel [63] it is not yet 
experimentally possible to perform them on selected pairs of atoms. However, 
proposals on utilizing marker atoms [46] or on quantum computing via global 
system control [67] exist and will allow to overcome this problem in the near 
future. 



17 



5. 1 Entanglement creation via coherent ground state collisions 



We exploit the two particle interaction term g J d 3 x^(x.)^^(-K)^(x.)^(-x.)/2 
from Eq. (4) whose origin lies in the atom-atom interaction potential for the 
controlled creation of entangled states. In optical lattices these nonlinear atom- 
atom interactions can be large [12], even for interactions between individual 
pairs of atoms. In addition we use a state selective optical lattice potential as 
introduced in Sec. 2 to provide control over the motional states of the atoms. 

Consider a situation where two atoms populating the internal states \a) and 
| b), respectively, are trapped in the ground states ipQ b of two potential wells 
V a ' b . Initially, at time t — —t , these wells are centered at positions x a and x b , 
sufficiently far apart (distance d = xi, — x a ) so that the particles do not interact 
(see Fig. 10a). The positions of the potentials are moved along trajectories 
x a (t) and x b (t) so that the wave packets of the atoms overlap for certain time, 
until finally they are restored to the initial position at t — r. This situation is 
described by the Hamiltonian 



/3=a,b 



I3\2 



2m 



+ V$(x a -x b ). (8) 



Here, V a > b (x a ' b - x a ' b (t)) describe the displaced trap potentials and is the 
interatomic interaction potential. Ideally, we want to implement the transfor- 
mation from before to after the collision, 

^(x a -x a )^ b (x b -x b ) -> e l ^(x a -x a )^ b (x b -x b ), (9) 

where each atom remains in the ground state of its trapping potential and 
preserves its internal state. The phase <fi = <p a + <p b + ab will contain a 
contribution ab from the interaction (collision) and (trivial) single particle 
kinematic phases 4> a and <p b . The transformation (9) can be realized in the 
adiabatic limit, whereby we move the potentials slowly on the scale given by 
the trap frequency, so that the atoms remain in their ground state. In this 
case the collisional phase shift is given by ab = J^dtAEit), where AE(t) 
is the energy shift induced by the atom-atom interactions, which in the case 
of purely s-wave scattering is given by 

AE(t) = J da#o (x - x a (t)) | 2 |^o (x - x b (t)) | 2 . (10) 

In addition we assume that \AE(t)\ <^ uj t so that no sloshing motion is 
excited. 

This process of colliding atoms by hand can be used to create entangled states 
between the two atoms. By moving a state selective optical lattice potential 
the wave function of each atomic qubit splits up in space according to the 
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atom 1 atom 2 




Fig. 10. We collide one atom in internal state \a) (filled circle, potential indicated 
by solid curve) with a second atom in state \b) (open circle, potential indicated by 
dashed curve)). In the collision the wave function accumulates a phase according to 
Eq. (9). a) Configurations at times t = ±r and b) at time t. 

internal superposition of states \a) and \b) as shown in cf. Fig. 10b. Only the 
wave function of the left atom in state \a) reaches the position of the second 
atom in state \b) and they will interact with each other. However, any other 
combination of internal states will not interact and therefore such a collision 
in a state selective optical lattice is conditional on the internal state. 

A maximally entangled Bell state can be created by carrying out a state 
selective collision with a phase ab = 7r and an initial superposition |+) oc 
\a) + \b) of each atom. After the collision this state is transformed into \a) |+) + 
\b) |— ) with |— ) oc \a) — \b) , i.e. |+) of the second atom is flipped to the 
orthogonal state |— ) conditional on the state of the first atom thus creating a 
Bell state. Identifying \a) = |0) and \b) = |1) we obtain for the logical states 
the mapping \i,j) — > (— 1)( 1_ *)J with i,j 6 {0, 1} (up to trivial phases) 

which realizes a two-qubit phase gate. This gate in combination with single- 
qubit rotations is universal i.e. they allow to implement any unitary dynamics 
on a quantum register [60]. 

5.2 Universal Quantum Simulators 

Building a general purpose quantum computer, which is able to run for exam- 
ple Shor's algorithm, requires quantum resources which will only be available 
in the long term future. Thus, it is important to identify nontrivial appli- 
cations for quantum computers with limited resources, which are available 
in the lab at present. Such an example is provided by Feynman's universal 
quantum simulator (UQS) [68]. A UQS is a controlled device that, operat- 
ing itself on the quantum level, efficiently reproduces the dynamics of any 
other many-particle system that evolves according to short range interactions. 
Consequently, a UQS could be used to efficiently simulate the dynamics of a 
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generic many-body system, and in this way function as a fundamental tool for 
research in many body physics, e.g. to simulate spin systems. 

According to Jane et al. [69] the very nature of the Hamiltonian available in 
quantum optical systems makes them best suited for simulating the evolu- 
tion of systems whose building blocks are also two-level atoms, and having a 
Hamiltonianifjv — J2 a + J2 a ^b that decomposes into one-qubit terms 
and two-qubit terms H^ ab \ A starting observation concerning the simula- 
tion of quantum dynamics is that if a Hamiltonian K = X^=i Kj decomposes 
into terms Kj acting in a small constant subspace, then by the Trotter for- 
mula e~ lKr = limm^oo {^r' lK ^ T l m e -' lK ^ T l m . . . e -iK S T/m\ we can a pp rox i m ate 
an evolution according to K by a series of short evolutions according to the 
pieces Kj. Therefore, we can simulate the evolution of an iV-qubit system ac- 
cording to the Hamiltonian Hn by composing short one-qubit and two-qubit 
evolutions generated, respectively, by and H^ ah \ In quantum optics an 
evolution according to one-qubit Hamiltonians can be obtained directly 
by properly shining a laser beam on the atoms or ions that host the qubits. 
Instead, two-qubit Hamiltonians are achieved by processing some given in- 
teraction Hq (see the example below) that is externally enforced in the 
following way. Let us consider two of the N qubits, that we denote by a and 
b. By alternating evolutions according to some available, switchable two qubit 
interaction Hq for some time with local unitary transformations, one can 
achieve an evolution 

n n n 

U(t = £ tj) = II Vj eM-iHt\)V] = n expH^r^) 
j=\ j=\ j=\ 

where t = X^=i tj, Vj = uf*' *®vf^ with Uj and Vj being one-qubit unitaries. For 

a small time interval U(t) ~ 1 - itY,]=iPjV j H^ ab) Vj + 0{t 2 ) with pj = tj/t, 

so that by concatenating several short gates U(t), U(t) = exp(-iH^'h) + 
0(t 2 ),we can simulate the Hamiltonian 

n 

for larger times. Note that the systems can be classified according to the 
availability of homogeneous manipulation, Uj = Vj, or the availability of local 
individual addressing of the qubits, Uj ^ Vj. 

Cold atoms in optical lattices provide an example where single atoms can 
be loaded with high fidelity into each lattice site, and where cold controlled 
collisions provide a way of entangling these atoms in a highly parallel way. 
This assumes that atoms have two internal (ground) states |0) = | [) and 
|1) = | |) representing a qubit, and that we have two spin- dependent lattices, 
one trapping the |0) state, and the second supporting the |1). An interaction 
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between adjacent qubits is achieved by displacing one of the lattices with 
respect to the other as discussed in the previous subsection. In this way the 
|0) component of the atom a approaches in space the |1) component of atom 
a + 1, and these collide in a controlled way. Then the two components of each 
atom are brought back together. This provides an example of implementing 
an Ising J2 a ^b = J2a ^ ® o"i a+1 - ) interaction between the qubits, where 
the a^'s denote Pauli matrices. By a sufficiently large, relative displacement 
of the two lattices, also interactions between more distant qubits could be 
achieved. A local unitary transformation can be enforced by shining a laser 
on the atoms, inducing an arbitrary rotation between |0) and |1). On the time 
scale of the collisions requiring a displacement of the lattice (the entanglement 
operation) these local operations can be assumed instantaneous. In the optical 
lattice example it is difficult to achieve an individual addressing of the qubits. 
Such an addressing would be available in an ion trap array as discussed in 
Ref. [70]. These operations provide us with the building blocks to obtain an 
effective Hamiltonian evolution by time averaging as outlined above. 

As an example let us consider the ferromagnetic [antiferromagnetic] Heisen- 
berg Hamiltonianif = J J2j= x , y ,z a j ® a j where J > [ J < 0] . An evolution 
can be simulated by short gates with Hq 1 ^ = 7<r 2 <g> a z alternated with local 
unitary operations 
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without local addressing, as provided by the standard optical lattice setup. 
The ability to perform independent operations on each of the qubits would 
translate into the possibility to simulate any bipartite Hamiltonians. 

An interesting aspect is the possibility to simulate effectively different lattice 
configurations: for example, in a 2D pattern a system with nearest neighbor 
interactions in a triangular configuration can be obtained from a rectangular 
array configuration. This is achieved by making the subsystems in the rectan- 
gular array interact not only with their nearest neighbor but also with two of 
their next-to-nearest neighbors in the same diagonal (see Fig. 11). 

One of the first and most interesting applications of quantum simulations is 
the study of quantum phase transitions [71]. In this case one would obtain the 
ground state of a system, adiabatically connecting ground states of systems in 
different regimes of coupling parameters, allowing to determine its properties. 
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triangular lattice 



square lattice 



Fig. 11. Illustration how triangular configurations of atoms with nearest neighbor 
interactions may be simulated in a rectangular lattice. 



6 Impurities and Atomic Quantum Dots 



A focused laser beam superimposed to a trap holding a quantum degenerate 
bose or fermi gas can form a tight local trapping potential, which can play the 
role of an impurity site, or an atomic quantum dot, coupled to a cold atom 
reservoir [72,73]. This opens the possibility of studying impurity physics or 
the analogue of quantum dot physics with quantum degenerate gases with 
controllable external parameters, and has several interesting applications, e.g. 
extraction of single atoms from an atomic BEC. Impurity potentials can be 
created as a local deformation of the original trap holding the atomic gas with 
tunnel coupling between reservoir and the dimple, or with a spin-dependent 
potential. In the latter case, atoms in the internal state \a) see only the large 
trap representing the reservoir, while atoms in state \b) are held by the lo- 
cal tight trap. Coupling of atoms occupying the atomic quantum dot to the 
reservoir of 6-atoms is achieved via a Raman transition (Fig. 12). Depending 
on the choice of the trapping potential for the reservoir atoms, we can realize 
configurations, where the dot is coupled to a reservoir of a ID (Luttinger), 2D 
or 3D quantum gas. Atoms occupying the dot and the reservoir will interact 
via collisional interactions. In particular, in the limit of a tight trap represent- 
ing the dot and strong repulsion between the atoms we can have a situation 
of a collisional blockade regime [73], where only a single atom can occupy the 
dot, in analogy to the Coulomb blockade in electronic quantum dots. In the 
following we will consider two specific examples to illustrate these configu- 
rations. Again our emphasis will be on the derivation of the corresponding 
Hamiltonians, while we refer to the original papers for a discussion of the 
system dynamics. 
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Fig. 12. Schematic setup of an atomic quantum dot coupled to a superfluid atomic 
reservoir. The Bose-liquid of atoms in state \a) is confined in a shallow trap V a (x). 
The atom in state \b) is localized in tightly confining potential V&(x). Atoms in state 
| a) and \b) are coupled via a Raman transition with effective Rabi frequency £1. A 
large onsite interaction Ubb > allows only a single atom in the dot. 

6.1 Atomic quantum dot coupled to superfluid reservoir: quantum engineering 
of a spin-boson model 



As our first example, we consider the configuration outlined in Fig. 12 where 
a superfluid reservoir of Bose atoms in internal state | a) is coupled via Raman 
transitions to the lowest vibrational state of an atomic quantum dot with 
atoms in state \b) [73]. The resulting dynamics will be described by a spin- 
Boson model [74], where the occupation of the dot (represented by a pseudo- 
spin) is coupled to the density and phase fluctuations of the superfluid. The 
distinguishing features of the model are: (i) the control of the spin-phonon bath 
coupling by the collisional parameters, and (ii) the realization of an ohmic or 
superohmic reservoir. As a result, the rich dynamics for the spin-Boson model, 
for example dissipative phase transitions, should be observable in the quantum 
dot dynamics. 

The collisional interaction of atoms in the two internal levels are described by 
a set of parameters g a/ 3 = Ama a plm with scattering lengths a a p (a, (3 = a,b). 
Following steps analogous to the derivation of the BHM Hamiltonian in Sec. 3 
we arrive at the Hamiltonian for the dot Hf, and its coupling to the reservoir 

H a b, 

H b + H ab = (-5 + g ab J dx |^(x)| 2 p a (x)) b r b 

+ Ej±bWb + J dxfi(* a (x)^(x)St + h.C.), (11) 

Here b is the bosonic destruction operator for 6-atom in the lowest vibra- 
tional state of the dot with wave function ^(x), while ^/ a (x) is the anni- 
hilation operator for a reservoir atom a at point x. The first term in (11) 
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contains the Raman detuning S and the collisional shift of the lowest vibra- 
tional state of the dot due to interaction with the reservoir with atomic density 
p a (x) = ^(x)^ f a (x). The second term is the onsite repulsion C/ 66 of the atoms 
occupying the quantum dot, and the third term is the laser induced coupling 
transferring atoms from the dot to the reservoir with effective two-photon 
Rabi frequency Q. 

At sufficiently low temperatures the reservoir atoms in \a) form a superfluid 
Bose liquid with an equilibrium liquid density p a . The only available excita- 
tions at low energies are then phonons with linear dispersion oj = v s \q\ and 
sound velocity v s . We write the atomic density operator in terms of the den- 
sity fluctuation operator IT, p a (x) = p a + II(x), which is canonically conjugate 
to the superfluid phase 0(x) with ty a ~ p a 1//2 e - ^. In the long wavelength 
approximation the dynamics of the superfluid is described by a (quantum) 
hydrodynamic Hamiltonian, H a = i> s ^ q Iql^q^qj which has the form of a col- 
lection of harmonic sound modes with annihilation operators of phonons 6 q . 
Since the excitations of a weakly interacting Bose-liquid are phonon-like only 
for wavelengths larger than the healing length £ > l b , the summation over the 
phonon modes is cutoff at a frequency u c = v s /t; ~ g aa Pa- 

In the collisional blockade limit of large onsite interaction U bb only states with 
occupation n b — and 1 in the dot participate in the dynamics, while higher 
occupations are suppressed by the large collisional shift. This situation and 
its description is analogous to the Mott insulator limit in optical lattices. 
Thus, the quantum state of a dot is described by a pseudo-spin-1/2, with 
the spin-up or spin-down state corresponding to occupation by a single or no 
atom in the dot. Using standard Pauli matrix notation, the dot occupation 
operator ¥b is then replaced by (1 + a z )/2 while ¥ — > a + . Furthermore, for 
the long wavelength wavevectors |q|Z& <C 1 with l b the dimension of the dot, 
the phonon field operators may be replaced by their values at x = 0. This 
gives the Hamiltonian 

H b + H ab = (- 5 - + ^»n(0)) a z + |(<7 +e -*(°> + h.c). (12) 

Here A ~ ^n]/ 2 is an effective Rabi frequency with n a the number of bosons 
inside the quantum dot wavefunction ipb- An important feature of this Hamil- 
tonian is that phonons couple to the dot dynamics via (i) the collisional 
interactions coupling to the density fluctuations 11(0), and (ii) to the su- 
perfluid phase 0(0) via the laser couplings. A unitary transformation H = 
S' _1 (i? a + Hb + H a b)S with S = exp(— a z i(f>(0)/2) transforms the phase fluctu- 
ations in the Rabi term to a frequency fluctuation which adds this contribution 
to the first term in (12). Thus the dynamics of the quantum dot coupled to 
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the superfluid reservoir gives rise to a spin-Boson Hamiltonian [74] 

# = E "q&i&q + + E \(&q + y " J0- X (13) 

with amplitudes of the total phonon coupling A q = \mqVg /2Vp a \ 1 ^ 2 (gabPa/rnv 2 — 1), 
as a sum of the collisional interaction and the transformed phase fluctua- 
tions. We note the following two features of the model. 

Control of spin-bath coupling by quantum interference: For a repulsive inter- 
species interaction, g a ^ > 0, both contributions in X q interfere destructively: 
the effect of the phonon excitation in a laser-driven transition a <-> b can be 
precisely cancelled by the change in the direct "elastic" interaction between 
the liquid and the b atoms, as described by the first term in Eq. (12). In a 
weakly interacting gas mv 2 s = p a g aa and thus the coupling constants A q vanish 
at g a b = gaa- At this special point, we have thus decoupled the dissipative 
phonon reservoir from the spin, i.e. we have stable superposition of occupation 
and non-occupation of the dot, exhibiting undamped Rabi oscillations of the 
quantum dot's occupancy. 

Ohmic and superohmic bath couplings: The system is characterized by the ef- 
fective density of states J{uS) = J2q K^i^— u; q ) = 2au s where a ~ (g a bp a /mVg — 
l) 2 is the dissipation strength due to the spin-phonon coupling and D = s the 
dimension of the superfluid reservoir. In the standard terminology [74], s = 1 
and s > 1 correspond to the ohmic and superohmic cases, respectively. 

For a discussion of the dynamics and the validity of the model, and specific 
atomic realizations we refer to [73] and [74]. 

6.2 A "Single Atom Transistor" with a ID optical lattice 

As a second example, we consider a spin-1/2 atomic impurity which is used 
to switch the transport of either a ID Bose-Einstein Condensate (BEC) or a 
ID degenerate Fermi gas initially situated to one side of the impurity [75]. 
In one spin state the impurity is transparent to the probe atoms, while in 
the other it acts as single atom mirror prohibiting transport. Observation of 
the atomic current passing the impurity can then be used as a quantum non- 
demolition measurement of its internal state, which can be seen to encode 
a qubit, \ip q ) = a\ |) + j3\ [). If a macroscopic number of atoms pass the 
impurity, then the system will be in a macroscopic superposition, [^(t)) = 
a|t)|0l(t)) + P\[)\4>i{t)) , which can form the basis for a single shot readout 
of the qubit spin. Here, \4> a (t)) denotes the state of the probe atoms after 
evolution to time t, given that the qubit is in state a (Fig. 13a). In view of the 
analogy between state amplification via this type of blocking mechanism and 
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Fig. 13. (a) A spin 1/2 impurity used as a switch: in one spin state it is transparent to 
the probe atoms, but in the other it acts as a single atom mirror, (b) Implementation 
of the SAT as a separately trapped impurity q with probe atoms b in an optical 
lattice. 

readout with single electron transistors (SET) used in solid state systems, we 
refer to this setup as a Single Atom Transistor (SAT). 

We consider the implementation of a SAT using cold atoms in ID optical 
lattices: probe atoms in state \b) are loaded in the lattice to the left of a 
site containing the impurity atom \q) , which is trapped by a separate (e.g., 
spin-dependent) potential (Fig. 13b). The passage of \b) atoms past the impu- 
rity q is then governed by the spin-dependent effective collisional interaction 
Hint = J2a U c s tU blboqlq a . By making use of a quantum interference mechanism, 
we engineer complete blocking (effectively U c s — > oo) for one spin state and 
complete transmission (U e s 0) for the other. 

The quantum interference mechanism needed to engineer U e s can be produced 
using an optical or magnetic Feshbach resonance [42,45,46], and we use the 
present example to illustrate Hamiltonians for impurity interactions involving 
Feshbach resonances and molecular interactions. For the optical case a Raman 
laser drives a transition on the impurity site, 0, from the atomic state &og£ |vac) 
via an off-resonant excited molecular state to a bound molecular state back 
in the lowest electronic manifold mj. |vac) (Fig. 14a). We denote the effective 
two-photon Rabi frequency and detuning by Q a and A CT respectively. The 
Hamiltonian for our system is then given by H — i/& + Hq, with 



H b = -JYfth + ^»E% - 1) ( 14 ) 

(ij) 3 

H = (™i<Zo-&o + h.cj - A a mlm a (15) 

cr 

Here H^ is a familiar Hubbard Hamiltonian for atoms in state \b); H describes 
the additional dynamics due to the impurity on site 0, where atoms in state 
| b) and \q) are converted to a molecular state with effective Rabi frequency 
Q a and detuning A a , and the last two terms describe background interactions, 
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U a/ 3 t a for two particles a, (5 e {g<7, which are typically weak. This model 

is valid for U a p, J, f2, A <C ujt (where oot is the energy separation between 
Bloch bands). Because the dynamics for the two spin channels q a can be 
treated independently, in the following we will consider a single spin channel, 
and drop the subscript a. 

To understand the qualitative physics behind the above Hamiltonian, let us 
consider the molecular couplings and associated effective interactions between 
the | q) and \b) atoms for (i) off- resonant (Q <C |A|) and (ii) resonant (A = 0) 
laser driving. In the first case the effective interaction between \b) and \q) 
atoms is U e s = U q b + Q 2 /A, where the second term is an AC Stark shift 
which plays the role of the resonant enhancement of the collisional interactions 
between \b) and \q) atoms due to the optical Feshbach resonance. For resonant 
driving (A = 0) the physical mechanism changes. On the impurity site, laser 
driving mixes the states &o^|vac) an d rn^\va,c), forming two dressed states 
with energies e± = (U qb )/2 ± (C/ ? 2 b /4 + tt 2 ) 1/2 (Fig. 14b, II). Thus we have two 
interfering quantum paths via the two dressed states for the transport of \b) 
atoms past the impurity. In the simple case of weak tunneling Q ^> J and U q b = 
second order perturbation theory gives for the effective tunnelling J eS = 
~ — j^q — ¥ (|e| fi) which shows destructive quantum interference, 
analogous to the interference effect underlying Electromagnetically Induced 
Transparency (EIT) [76], and is equivalent to having an effective interaction 
U cS — > oo. 

In Ref. [75] the exact dynamics for scattering of a single \b) atom from the 
impurity is solved exactly, confirming the above qualitative picture of EIT- 
type quantum interference. Furthermore, in this reference a detailed study of 
the time-dependent many body dynamics based on the ID Hamiltonian (16) 
is presented for interacting many-particle systems including a ID Tonks gas. 



7 Conclusions 



The cold atom Hubbard toolbox presented in this paper opens a range of novel 
and exciting prospects for controlling and manipulating strongly correlated 
atomic systems. The resulting standard and exotic many body Hamiltonians 
can be used to simulate otherwise untractable systems in condensed matter 
physics. The quantum simulation of such systems becomes feasible with cold 
atoms and provides a valuable tool which can be used in addition to standard 
analytical and numerical methods for analyzing many body systems. In addi- 
tion we have shown that they also pose new challenges like e.g. understanding 
and utilizing the unitary dynamics that can be realized in optical lattices. 
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Fig. 14. (a) The optical Feshbach setup couples the atomic state &Q<fJ-|vac) (in a 
particular motional state quantized by the trap) to a molecular bound state of 
the Born-Oppenheimer potential, mi- 1 vac), with effective Rabi frequency Q a and 
detuning A CT . (b) A single atom passes the impurity (I— ►III) via the two dressed 
states (II), |+) = P qt\vac) + m^vac) and |— ) = frj^lvac) — mj-|vac). Quantum 
interference between the paths gives rise to an effective tunnelling rate J c s,a- 
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